Generating-function approach for bond percolations in hierarchical networks 
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We study bond percolations on hierarchical scale-free networks with the open bond probability 
of the shortcuts p and that of the ordinary bonds p. The system has a critical phase in which 
the percolating probability P takes an intermediate value < P < 1. Using generating function 
approach, we calculate the fractal exponent ip of the root clusters to show that ijj varies continuously 
with p in the critical phase. We confirm numerically that the distribution n s of cluster size s in the 
critical phase obeys a power law n s oc s~ T , where r satisfies the scaling relation r = 1 + ip -1 . In 
addition the critical exponent f3(p) of the order parameter varies as p, from f3 ~ 0.164694 at p — 
to infinity at p = p c = 5/32. 
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I. INTRODUCTION 

Dynamics on and of complex networks have been one 
of the focuses of attentions since late 1990s 0-0]. Real 
networks, e.g., WWW, Internet, food- web, often have 
complex properties such as scale-free degree distribution 
[H, small- world property [f|, etc., to demand more ex- 
tensive framework of statistical physics to investigate the 
interplay between dynamics and such network topology. 
Many analytical and numerical works about the effects 
of network topology on processes such as percolation, in- 
teracting spin systems, epidemic processes, have been re- 
ported jg]. 

Among these issues, unusual phase transitions of per- 
colations and spin systems on some networks have at- 
tracted our current interests For example, the 
percolations on some growing network models undergo 
an infinite order transition with a Berezinskii-Kosterlitz- 
Thouless (BKT)-like singularity: (i) the relative size of 
the largest component vanishes in an essentially singular 
way at the transition point, so that the transition is of 
infinite order, and (ii) the mean number n s of clusters 
with size s per node (or the cluster size distribution in 
short) decays in a power-law fashion with s, 



n s oc s 



(1) 



in a finite region below the transition point where no gi- 
ant component exists @~QjI • A similar non-ordered phase 
with some power-law behavior, a critical phase, has also 
been observed in bond percolations on the enhanced bi- 
nary tree [l7rjT9| . which is one of nonamenable graphs 



(NAGs) [2fj,[21(. The system on a NAG takes three dis- 
tinct phases according to the open bond probability p as 
follows: (i) the non-percolating phase (0 < p < p cl ) in 
which only finite size clusters exist, (ii) the critical phase 
(Pel < P < Pc2) in which there arc infinitely many infinite 
clusters, and (iii) the percolating phase (p C 2 < p < 1) in 
which the system has a unique infinite cluster. Here in- 
finite cluster means a cluster whose mass diverges with 
system size N as with < <f> < 1. To profile the crit- 
ical phase it is useful to calculate the fractal exponent 
ip-L defined as s max cx N^ h , where s max is the mean size 
of the largest components in the system with N nodes. 
Note that corresponds to df/d for percolating clus- 
ters having the fractal dimension df on ti-dimensional 
Euclidean lattices. Recent paper [17| has shown nu- 
merically that the above phases are characterized as (i) 
V'lO?) = for p < pd, (ii) continuously increasing of 
iphip) (0 < iph(p) < 1) with p, where n s also behaves as 
([T]) with p-dependent r satisfying 



(2) 
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for pd < p < p c2 , and (iii) iph(p) = 1 for p > p C 2- 
The scaling relation ([2]) indicates that Vl plays a role 
of the natural cut-off exponent of n s as shown in the 
growing random tree in which p c i = and p C 2 = 1 pd| . 
In general the growing random networks are considered 
to have p c \ = with finite p C 2 fMwty . 

There exist other systems having a similar phase. They 
are in a special class of hierarchical scale-free networks, 
called (decorated) (u, u)-flower introduced comprehen- 
sively in [22|, [23|. Berker et al. [l6[ have studied bond 
percolations on the decorated (2,2)-flower by renormal- 
ization group (RG) to show the existence of a critical 
phase (as known as the partially ordered phase 24|). 
where RG flow converges onto the line of nontrivial stable 
fixed points. But we have little knowledge about physical 
properties of the critical phase, i.e., how it is critical. 
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FIG. 2: (Color online) Phase diagram of the bond percolation 
on the decorated (2,2)-flower [l6l |. The blue-solid and the red- 
dashed lines denote p*(p) (stable fixed point) and p c (p) (un- 
stable fixed point), respectively. The shaded region represents 
the critical phase. 



FIG. 1: (Color online) Construction of the (2, 2)-flower F n 
and the decorated (2,2)-flower F n . (a) Each bond is replaced 
by two parallel paths consisting of two bonds each at next 
generation, (b) The flower F n +i of the n + 1-th generation is 
obtained by joining four copies of F n . (c) Realization of F n 
with n = 0, 1, 2, 3. (d) Realization of F n with n = 0, 1, 2, 3. 
Fn is obtained by adding the shortcuts (orange-dashed line in 
bond replacement (a)) to F n . The shortcuts are not replaced 
by others in each iteration. 



In this paper, we investigate bond percolations on the 
decorated (2, 2)-fk>wer with two different probabilities p 
and p, which are the open bond probability of the short- 
cuts and that of the ordinary bonds, respectively. Here 
we adopt a generating function approach to calculate the 
fractal exponent, the cluster size distribution, and the 
order parameter for an arbitrary combination of p and 
p. to reveal a complete picture about the phases of this 
model. Our calculations show (i) the fractal exponent 
ip and ft of the order parameter depend on the existing 
probability p of the shortcuts, and (ii) n s is a power-law 
at all the point in the critical phase, and its exponent r 
also depends on p. 

The organization of this paper is as follows: In SccllTl 
we introduce the (2,2)-flower and the decorated (2,2)- 
flower, and briefly review the pre vious studies for the 
percolation on the flowers (lH |22|. [23j . In Sec lIIIl we 
introduce the generating functions, and derive those re- 
cursion relations to calculate the order parameter, fractal 
exponent, and cluster size distribution. The main results 
are presented in Sec lIVl and Section [V] is devoted to sum- 
mary. 



II. MODEL 

In this section we briefly introduce the (2,2)-flower and 
the decorated (2,2)-flower. The (2,2)-flower F n of the 
n-th generation is constructed recursively as illustrated 
in Fig[TJ At n = 0, the flower Fq consists of two nodes 
connected by a bond. Hereafter we call these nodes roots. 
For n > 1, F n is obtained from F n _i, such that each 
existing bond in F„_i is replaced by two parallel paths 
consisting of two bonds each. The decorated (2,2)-flower 
F n , a variant of the (2,2)-flower F n , is given by adding 
shortcuts to F n , as illustrated in FigQJd). 

The network properties of these two flowers have been 
reported in [22|, [23| ■ The number of nodes N n of the n-th 
generation is N n = 2(4" + 2)/3 and the degree distribu- 
tion has a scale-free form P(k) oc fc -7 with 7 = 3 for both 
flowers. One of the important differences between F n and 
F n appears in those dimensionality. Since the diameter 
L n of F n is L n = 2", the dimension d of the underlying 
network defined as N n oc is 2. On the other hand F n 
is known to have small- world property L n ~ ln N n corre- 
sponding to d ~ oo. In addition F n has a high clustering 
coefficient C ~ 0.820, in contrast to C = for F n . 

In the present work we consider the bond percolation 
on F n with the open bond probability p of the bonds con- 
stituting F n (the ordinary bonds) and that of the short- 
cuts p being given independently. The standard bond 
percolation on F n and F n are given by setting p = and 
p = p, respectively. Note that the latter is also given by 
setting p = because F n with p = p and with p = 
is exactly the same. 

The phase diagram is solved exactly by RG technique 
[3, H2|. Let P (n) be a probability that both roots are 
in the same cluster of a bond configuration on F n with 
fixed p. The initial value is set to — p. In the 

large size limit, the system is regarded as percolating if 
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FIG. 3: Schematic for generating functions (a) T n (x), (b) 
S„(x,y), and (c) U n (x). The open circles represent the root 
nodes. 



P := linin^oo is nonzero. In this sense P or pW 
is called the percolation probability. Since p( n ) is given 
recursively as 



P (n+l) =i_(i_^ (l-(p(")) 2 y 



(3) 



one obtains the flow diagram from the solution of the 
equation (Fig|2~]). For < p < p c — 5/32, there are two 
nonzero stable fixed points, P — p*(p) < 1 and P = 1, 
corresponding to the partially ordered phase and the or- 
dered phase, respectively, and one unstable fixed point 
between the two giving the phase boundary, P = p c (p)- 
For p > p c , on the other hand, there is only one sta- 
ble fixed point at P = 1, so that the system is always 
percolating. 

Two special cases, p = and p = p, were investigated 
in [23]. Here let us recall their results briefly. For the 
case of p — 0, i.e., the standard bond percolation on F n , 
Ecj-® gives the critical point p c (p — 0) = (V5 — l)/2. A 
simple RG argument then gives the critical exponents at 
p c ; the exponent j3 ~ 0.164694 of the order parameter, 
i.e., the fraction of the largest component Poo ~ (p—pc) 13 , 
and v ~ 1.63528 of the correlation length £ ~ \p — p c \~ v } 
which are close to those of the two dimensional regular 
systems. On the other hand the same argument for p = p 
gives the infinite order transition, i.e., /3 — > oo. 

According to the above definition of percolation the 
percolating cluster should include both of the root nodes. 
It is then convenient to consider the mean size (so) n of 
the cluster including both roots (referred to as the root 
cluster) on F n instead of s max (A r „) to characterize the 
criticality: 



(4) 



where tp is the fractal exponent for the root cluster. Note 
that ijj behaves essentially the same as tJjl for the growing 
random trees UM. 



III. GENERATING FUNCTIONS 

In this section we describe how to utilize generating 
functions to calculate the fractal exponent ip and the 




FIG. 4: Possible diagrams contributing to (a) T n +i(x), (b) 
S n +i(x, y), (c) U n +i(x). The roots (open circles) are not 
counted in the generating functions, so that nodes connecting 
two F n s (closed circles) are taken into account by multiply- 
ing x or y. For example the first diagram of (a) represents 
x 2 T^(x) and the fourth diagram of (b) xyT n (x)T n (y)Sn(x, y). 



cluster size distribution n s (p) with p fixed. First let 
us consider the bond percolation on F n with open bond 
probability p with p = 0. We introduce three basic quan- 
tities on F n : the probability tC (p) that both roots arc 
connected to the same cluster of size k, the probability 
s£j (p) that the left (right) root is connected to a cluster 
of size k (I) but these clusters are not the same, and the 
mean number (p) of clusters of size k to which neither 
of the roots is connected. For the sake of convenience the 
roots arc not counted in the cluster size k or I for ijj. (p) 



defined as 



and Sfc"/ (p) . The corresponding generating functions are 



T n {x) 
S„(x,y) 

U n {x) 



E4 n) (p)* fe , 



k=0 
oo oo 



k=0 1=0 
oo 



(5a) 
(5b) 
(5c) 



k=0 
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The self-similar structure of F n allows us to obtain the 
recursion relations for the above generating functions 



T n+1 (x) = x 2 T*(x)+4x 2 T*(x)S n (x,x) 
S n +i(x,y) 



(6a) 



U n+ i(x) 



+2xT 2 (x)S 2 n (x,l), 
S 2 (x,l)S 2 (l,y) 
+2S n (x, y)S n (x, l)S n (l, y) 

x [xT n (x) + yT n (y)} 
+ 2x ySn { x , v) T n {x)T n (y) 
+Sl{x,y) [x 2 T 2 (x)+y 2 T 2 (y)]A^) 
W n {x) + 2xS 2 n {x,l) 1 (6c) 



as illustrated in FigJ3] and FigJU The initial conditions 
are given as Tq(x) — p, So(x,y) = q = l—p and Uq(x) = 
0. 

It is convenient to rewrite these recursion formulas in 
terms of functions of single variable x only. In order to 
this we introduce new functions V n (x) = S n (x,x) and 
R n (x) = S n (x, 1) to obtain 



T n+1 (x) = &[T n (x),V n (x),R n (x),x] 

= x 2 T*(x) + 4x 2 T*(x)V n (x) 
+2xT 2 (x)R 2 n (x), 
V n+1 (x) = f [T n (x),V n (x),R n (x),x] 

= RUx)+AxT n (x)V n {x)R 2 n (x) 
+4x 2 T 2 (x)V 2 (x), 
R n+1 (x) = >M[T n (x),R n (x),x} 

= Rl(x)[l + xT n (x)] 2 , 
U n+1 (x) = ^[R n (x),U n (x),x] 

= 4U n (x)+2xR 2 n {x). 



(7a) 

(7b) 
(7c) 
(7d) 



Indeed it is this form that enables us to obtain the solu- 
tions for large n numerically. 

Now the construction of the recursion relations for 
F n with p fixed is straightforward. Let T n (x), V n {x), 
R n {x) and U n (x) denote the corresponding generating 
functions on F n . By using the above formula ([7aj) - (|7d[) 
with these functions one can construct the generating 
functions T n +i(x), V n +i(x), R n +i(x) and U n +i(x) on the 
decorated (2,2)-fiower of the next generation without the 
shortcut directly connecting the roots: 

T n+ i{x) = ^\ i f n (x),V n (x),R n (x),x\, (8a) 

V n+1 (x) = y[f n (x),V n (x),Rn(x),x], (8b) 

R n+1 (x) = <%[f n {x) 7 R n (x) 7 x] 7 (8c) 

U n+1 (x) = w\R n (x),U n (x),x] . (8d) 



The flower F n+ i is made by adding the shortcut to the 
intermediate one with probability p and one thus obtains 



f n+1 (x) = T n+1 (x) +pV n+ i(x), (9a) 

V n+1 (x) = qVn+x(x), (9b) 

R n+ i(x) = q~R n+ i(x), (9c) 

U n+1 {x) = U n+1 {x), (9d) 



where q = l—p. The initial conditions are given as 
T (x) = p, V~o( x ) = Ro(x) = q and Uo(x) = 0. One can 
easily check the probability conservation T n (l) + V n (l) = 
1 by the iteration (J9j) . 

Once these generating functions are obtained one can 
evaluate various quantities of the present interest. For 
example the percolation probability P'"' is given as 



P [n) =f n (l), (10a) 

Q („) = 1 _ p(n) = 1} = = (1Qb) 



Note that (j3|) is re-obtained by putting a; = 1 to fl§J) 
and using ([TU)) . The mean number of the root cluster 

(or the order parameter P^\p)) and the cluster 



size distribution n j 1 (p) on F„ are given as 



(so)» 
P^Hp) 



7^(1) + ^(1), 

(so)n _ 



~(n) 



GO 



(11) 
(12) 

(13) 



where the prime denotes the first derivative with respect 
to x, and we put r„ = T^(l)/N n and a n = V^(l)/N n . 

It is useful to consider the derivatives of recursion 
relations (|§al) - (f9"c)) for evaluating P<£/{p). By noticing 
V^(l) = 2_R^(1) we obtain the recursion relations for r„ 
and fT„ as 
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o~ n +l 

T~n+1 



4g(Q(")) 2 (l + p(")) 



A/„ +1 V2(l + P (n) ) [(-P ( " ) ) 2 +pQ ( " ) (l + ^ ( " ) )] 4[l-g(Q(")) 2 (l + p(™))] 

i f 4gp(™)(g(™)) 2 (i + p(™)) 



(for n > 1), 



N n+1 \2PW [2 - PW - 2g(Q(")) 2 (l + pW)] 

±gQ(l + P) 2 qQ 2 (l + P) A /o 

|(1 + P) [P 2 +pQ(l + P)] l-gQ 2 (l + P 



(14) 
(15) 



where we recall P = lim„_ i . 00 P^ and Q = lirrin^oo Q^ n \ 
Note that this expression is an extension of Eq.(31) in 

M. 



IV. RESULTS 

To profile the critical phase we calculate the fractal 
exponent if. In the ordered phase we have trivially if = 1. 
Otherwise the fixed points P(< 1) of the RG equation ^ 
satisfy q(l — P)(l + P) 2 = 1. The recursion relation (fTS"|) 
is then reduced to 



Cn+l 



; P 1 - a 



(16) 



where a = (1 — P)/(l + P). By using the largest eigen- 
value A(P) of the above matrix, 

1 



A(P) = - [(3 - 2a) + v/l-4a(l-2P) + 4a 2 J , (17) 

we can calculate the fractal exponent ip on the fixed 
points in the same way as [22j does: 



ip{P) = 1 + 



lnA(P) 
In 4 



(18) 



The p-dependence of ip is shown in Fig(5] We find that 
(i) for p < p c = 5/32, ip on the (un) stable fixed points in- 
creases (decreases) with increasing p, and (ii) for p > p c , 
ip is equal to one irrespective of both p and p, which 
means that the system is always in the percolating phase. 
Let us consider the p-dependence of if with p < p c fixed. 
In the critical phase (p < p c {P)) the percolation probabil- 
ity P( n > goes to p*(p) and the exponent ip = ip{p*{p)) is 
constant in this region. At the critical point (p = p c (p)) 
the fixed point is P = p c (p) itself and thus ip discontin- 
uously changes to ip(Pc(P)), and jumps again to one for 
the percolating phase (p > p c (j>)). This behavior can be 
also confirmed directly by evaluating the 7V n -dependcncc 
of (so)n numerically (not shown). This result indicates 
that the probability p of the shortcuts essentially deter- 
mines how the system is critical in the partially ordered 
critical phase. On the other hand, for the standard bond 
percolation on the decorated (2,2)-flowcr (p = p), the 
fractal exponent ip varies continuously with open bond 
probability p as observed on a NAG [17| . 




0.15\ 0.2 

5/32 



FIG. 5: (Color online) Fractal exponent ip on the stable fixed 
points p*(p) (blue-solid line) and unstable fixed points p c (p) 
(red-dashed line). 



We also evaluate numerically the recursion relations 
(f9aj) - (|9d[) to obtain the cluster size distribution n s on F n 
given by (|13|) . We should observe a power-law behavior 
of n s in the whole region of the critical phase. To check 
this behavior we assume a finite size scaling form 



n s (N) 



(19) 



where ip is the fractal exponent obtained above and the 
scaling function f(x) behaves as 

. ] rapidly decaying func. for x > 1, . . 

f( x ) ~ S . . f „ , I 20 ) 

I constant for x -C 1, 



and t satisfies the scaling relation |ll|, [lj 

T = l+i>-\ 



(21) 



As discussed in [FL( a plausible argument leads us to the 
relation as follows: By assuming n s oc s~ r in the critical 
phase, a natural cutoff s max of the cluster size distri- 
bution (a natural cutoff of the degree distribution was 
introduced in (25[) is given as 



JV 



n„ds ~ 1 — > 



cx AT' 



(22) 



Here we emphasize that s max plays just a role of char- 
acteristic cutoff of the distribution and does not need to 
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FIG. 6: (Color online) Finite size scaling for cluster size 
distribution n 3 at (a) stable fixed point (p = 0.130302) and 
(b) unstable fixed point (p — 0.517492), and (c) p = 0.3, 
for p = 0.1. Here ip and r in (c) are given by those at the 
stable fixed point. Insets show raw data of n s . We set (a) 
n = 9, fO, If, 12, 13, 14, (b) n = 5, 6, 7, 8, 9, f 0, and (c) n = 
8, 9, 10, 11, 12, 13, from left to right. 



be strictly the mean size of the largest clusters. In this 
sense we can replace the fractal exponent tph of s max with 
-0 of (sq), so that the relation (pH]) follows. Indeed the 
largest cluster not containing the roots is expected to the 
one containing the top-most (or, equivalently, bottom- 



a? 




FIG. 7: (Color online) Fraction of the largest components 
-Poo(p) on F n with p = 0, 0.1, 0.15, and 5/32(= p c ) (from right 
to left). Here n is taken to 10 6 (» 1). 
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FIG. 8: p-dependence of the critical exponent f3 on the phase 
boundary p = p c (p). Inset shows that of /3 _ . 



most) node in the right figure of FigJ^c) contributing 
the second term of (|6c)) . The mean size of the cluster is 
proportional to V^(l) and so to (so). Therefore one can 
conclude that the characteristic size grows with N not 
faster than N*. 

Our finite size scaling for n s is indeed well fitted on 
both stable and unstable fixed points as shown in Fig|51 
Note that for the sake of (j!?Tj) no fitting parameter re- 
mains. The scaling also works at any p in the critical 
phase (FiglHlJc)), but the convergence is not so rapid as 
on the fixed points. 

Finally we iterate Ea.lfFof numerically to obtain the 
order parameter P^\p) on F n . The result for n = 10 6 is 
shown in Fig(7] The initial growth of the order parameter 
becomes moderate with increasing p. To examine the 
critical exponent /3 on the phase boundary p = p c {p), we 
follow the scaling argument in [22j to obtain 



Pip) 



lnA(p c (p)) 
'lnA(p c (p))' 



(23) 
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where 



SUMMARY 



HP) 



gp(n+l) 



dP(r 



4<?P(1 - P 2 ) = 2(1 -a). (24) 



Figure [8] shows the p-dependence of /?. We find that 
j3 increases continuously with p, from f3 = 0.164694 at 
p = to j3 = oo at p = p c = 5/32. At p = p c we expand 
Eq.Q near p c (p c ) = 1/3 to obtain 



^p(n+i) ^ ^p(n) _|_ _(Ap(")) 2 , 



(25) 



where AP^™) = P^ n > —p c (p c ). We can estimate the solu- 
tion for small AP(°) = p — p c (p c ) > as 



AP W ~ AP (0) 



9n 



( A p(°)) 5 



(26) 



which is correct as long as the second term in the r.h.s. 
is much less than the first one, or equivalently, n is much 
less than n* ~ l/AP^ = |p - p c \~ x . For n > n*, pW 
goes to 1 rapidly and so we obtain 



P. 



A™ 



cxp 



const. 
P-Pc 



(27) 



where A = X(p c (p c )) given by Eq. ([T7|) . We thus find an 
essential singularity in the order parameter at p = p c . 

Note that j3 is apparently related to ip through X(p c (p)) 
as shown in [22[ (see Eqs. tfLS)) and (|23|> V It is, however, 
not the case in the (off-boundary) critical phase where 
the nontrivial stable fixed points p* (p) dominate the crit- 
icality while the order parameter vanishes there. 



We have investigated bond percolations on the deco- 
rated (2,2)-flower with two different probabilities p and 
p. Our generating function approach has revealed that 
the system is in the critical phase for p < p c {P) and 
p < p c = 5/32. We have evaluated the fractal exponent 
ip and confirmed the power-law behavior of n s in the crit- 
ical phase as well as those dependence on p and p and 
the validity of the scaling relation r = 1 + ip -1 . 

We have also examined the critical exponent /3 in the 
percolating phase and found that j3 also varies as p from 
/3 ~ 0.164694 at p = 0, where the network is two- 
dimensional-like, to (3 = oo at p = pc, where the di- 
mensionality of the underlying network is infinite. 

It is only at p = p c that the percolation on the deco- 
rated (2,2)-flower shows an infinite order transition with 
the BKT-like singularity as percolations on growing net- 
works do 0-E3] . The finiteness of /3 for p < p c suggests 
that the existence of some critical phase adjacent to the 
normal ordered phase is not enough for the network to 
have such an essential singularity in the order parameter 
and thus an infinite order phase transition. At present 
we have none of the key to reveal necessary conditions for 
the existence of the BKT-like singularity. Further study 
would be required to clarify the relation between these 
interesting properties of the phase transition. 
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